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A Comparison of the Fractal and JPEG Algorithms 
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A proprietary fractal image-compression algorithm and the Joint Photographic 
Experts Group (JPEG) industry standard algorithm for image compression are 
compared. In every case , the JPEG algorithm was superior to the fractal method 
at a given compression ratio according to a root-mean -square criterion and a pea 
sign a I- to- noise crite ri on . 


I. Introduction 

Fractal image compression has attracted much public- 
ity in recent years. It has been suggested that one can 
achieve compression ratios of the order of thousands to 
one by the application of fractal algorithms. Some re- 
searchers successfully generated certain images, with very 
small databases, by using fractal algorithms. These im- 
ages generally consisted of natural objects, and the mem- 
ory requirement for, e.g., a realistic looking tree, was about 
120 bytes. However, the applicability of these methods to 
general image compression and the achievement of the phe- 
nomenal compression ratios of thousands to one have been 
viewed with general skepticism. In order to make a com- 
parative study of the fractal versus Joint Photographic Ex- 
perts Group (JPEG) standard algorithms, the authors sent 
ten images to a vendor, Iterated Systems Inc., Norcross, 
Georgia. These images were compressed by their propri- 
etary fractal algorithms, and reconstructed using their de- 
compression package. The compression ratios were from 
about five to one to twenty to one. The mean square er- 
rors and peak signal-to-noise ratio (SNR) were compared 
to the corresponding ones for the JPEG algorithms at the 
same compression ratios. The latter approach proved to 
be superior in every case according to these criteria. 


The theoretical foundation and the practical implemen- 
tation of the fractal method for image generation is de- 
scribed in Section II. The relevant portion of the JPEG 
algorithm is briefly described in Section III, and finally in 
Section IV the results of the comparative study are given. 


II. Fractals 

There are various ways of defining fractals. The frame- 
work proposed by J. Hutchinson [1] has been the most 
successful approach for the study of fractals. To describe 
this method let S = {Si, ■ • • * S n ) be a finite set of affine 
transformations of IV. This means that if x G IV, then 
the effect of Si = (A iy Vi) on x is given by 


Si{z) = Mx) + Vi 

where Ai is a linear transformation of IV and € IV It 
will be assumed that Af s are nonsingular and Si s are con- 
tracting, i.e., || Si(x) - ( 3 /) || < \\ x - 2 /|| for all x,y G IV- 

Then the affine transformation Si x • • ■ Si k is also contract- 
ing and, therefore, has a unique fixed point that will be 
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denoted by Fi lf ... t i k . The fractal set S ) associated to S 

is, by definition, 

T(S) = closure({F lli ... >l -J all 0 < ii, • . • , i k < n 
and k = 1,2,3, • ■ •}) 

F(S) is a compact subset of R q . The self- similarity prop- 
erty of fractals is expressed by the fundamental equation 

i= 1 

In fact, 

Theorem 1. From [1], F(S) is the unique compact set 
K C R ? with the property 

n 

K = \JSi(K) ( 1 ) 

1=1 

Notice that each Si(K) is a replica of K , so that Eq. ( 1 ) 
does indeed express the self-similarity property of fractals. 
This characterization of fractals is also important in prac- 
tical applications. 

Example 1. Let K be a convex polygon in R 2 with ver- 
tices vi,* and let Sx,*",S n be the affine transfor- 
mations Si(x ) = -j- where 0 < a; < 1 . If 

the a,’s are not too small, then (JS,(K) = K and, con- 
sequently, F(S) = K by Theorem 1 of fractal sets. Thus, 
every convex polygon can be realized as an T{S) for some 
$. On the other hand, it is not hard to see that the bound- 
ary dK of the convex polygon K cannot be realized as a 
fractal set. 

To generate a fractal image, one starts with a set 
S = {Si ,"•,£«} of affine transformations. For every 
product S it •"Si 2 S il = S = ( A,v ) of length / < jfc, for 
some pre-assigned value k y one computes its unique fixed 
point F itr -- } F it = F = (7 - A)~ l (v). Note that since S 
is contracting, I- A is invertible. The coordinates F\ and 
F 2 of F are multiplied by a normalizing factor N, and then 
quantized to the nearest integer to obtain a point p with 
integral coordinates (pi,p 2 ) € Z 2 . A point q — (71,72) 
on the screen is white (black) as it is (or is not) of form 
(Pl»P 2 ), as described above. 

The fractal set F(S) i thus constructed, corresponds to 
a black-and-white image. In order to introduce grey levels 


into .F(«S), one would like to use the density of the points 
{F {l , • - , F ik ] as a measure of the brightness of the pixels. 
To do so, it is convenient to regard the procedure for the 
generation of a fractal set as a random process. In fact, let 
Pi» ” ■ .Pn he positive real numbers such that Hpi = 1, and 
x G R 9 . Consider the random process X where x — ► Si{x) 
with probability pi. This process has a unique stationary 
distribution p. The stationarity condition is expressed by 
the equation 

P = E piSi*(p) ( 2 ) 

where S^(p) is the transform of the measure p under the 
affine transformation S*. Note that Eq. ( 2 ) is a more 
precise version of the fundamental self-similarity prop- 
erty expressed by Theorem 1 . The measure p is the 
mathematical representation of a grey-scale image on the 
screen. The support of the measure p is the fractal set 
J-(S) y and is independent of the choice of positive num- 
bers pi. Furthermore, if £* denotes the discrete probabil- 
ity measure naturally assigned to the set {F tl , - - * , iqj / < 
k and all *1 , ■ ••,*/}, then 

Sk — ► p weakly 

if all pi = 1/n. 

By introducing an alternative method for generating 
the fractal set F(S) y one can take advantage of the prob- 
abilities { pi } in actual image generation. Consider a 
realization of the random process X. This can be in- 
terpreted as follows: Using a random-number generator, 
one generates a sequence of integers {*1 , where 

l < h < n and the integer j is chosen with probability 
Pj • A realization of the process X is then the sequence of 
points S{ l (z),Si 7 Si 1 (z) } - • Let Alt be the discrete prob- 
ability measure assigned to the set {S (l (z), S^S^fc), •• • , 
Sit-'SiMz)}. Then, one can show by standard argu- 
ments that 

Theorem 2 , With probability 1, Mk — ► p weakly. 

After possibly multiplying by a factor N and quantiz- 
ing, one can regard a point Sq • • • Si 2 Si l ( z ) in a realization 
of the process X as the point = 7 = (<71,72) € Z 2 , 

One can make a histogram of the number of times the 
points of the integer lattice Z 2 are generated in this fash- 
ion. Grey levels are accordingly assigned to the points so 
that the points generated more often have higher intensity 
(are whiter) than those generated fewer times. The addi- 
tional input consisting of the positive numbers {pi , - • - , p n } 
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allows one to have better control over the density of the 
points, i.e., the grey levels. 

In order for any procedure to have practical applica- 
tions, it must be stable relative to the variation of the 
parameters. In the case of fractals, one would like the fun- 
damental property of self-similarity, as expressed by The- 
orem 1, to have the necessary stability. This means that 
if the transformations Si are such that U»=i Si(K) is ap- 
proximately identical with I< } then the fractal set T(S) is 
also approximately identical with K. It is not difficult to 
establish, in a quantitative sense, the validity of this sta- 
bility, which yields a slight generalization of Theorem 1. 
It is this mild generalization of Hutchinson s theorem that 
has been widely publicized by M. Barnsley as the Collage 
Theorem [2]. 

Just as in Example 1, where Theorem 1 was used to 
generate convex polygons as fractal sets, other natural- 
looking objects were generated by appealing to the ap- 
proximate version and the experimental insight into the 
effects of the variation of the parameters on fi. In actual 
applications, it is also convenient to decompose the matrix 
Ai as a product 

/ 1 a\ (ba 0 \ ( cos 9 -sin 0 \ 

^ y 0 1 y \ 0 6'a J ^ sin 9 cos 9 ) 

where |66'| = 1 and 0 < a < 1. The advantage of using 
this decomposition is that the effects of the parameters 
can be more readily understood. The matrix 

f cos 9 - sin 9 \ 

y sin 9 cos 9 ) 

is a rotation through angle 9 . The matrix 

(ba 0 \ 

\0 b'a) 

represents scaling by factors ba and- b a , and 



can be regarded as a twist. The numbers 6, 6', and a should 
be chosen such that |6a|, |fc'a| < 1. By specializing the 
parameters 6, fc r , and a to 1, 1, and 0, respectively, one can 


already generate a number of very complex and natural- 
looking images. The interested reader should view these 
points as hints or general guidelines for experimentation 
with fractal image generation. 

Using the fractal method, a number of images were gen- 
erated by one of the authors and others (see, e.g., [2] and 
[3]). Fractal methods have been used recently at JPL for 
computer simulation of certain images. Since these im- 
ages were generated by storing only the parameters of a 
few affine transformations, one would like to believe that 
the fractal method can be used for efficient data, or more 
specifically, image compression. The systematization of 
the fractal method so that it becomes applicable to gen 
eral image compression has been attempted by a num- 
ber of researchers. One approach is to try to find a set 
5 _ {S ir --,S n } of affine transformations and positive 
numbers {p t } such that the corresponding stationary dis- 
tribution is a good approximation to a given image. While 
in theory this is possible, the number of the affine transfor- 
mations may be so large that the result will have no practi- 
cal value. Furthermore, since most objects do not exhibit, 
even remotely, the self-similarity property that is an es- 
sential feature of fractals, this direct approach is probably 
a futile one. By segmentation of an image, one will have 
better control over the choice of the affine transformations. 
However, the startling compression ratios of thousands to 
one will not be achievable. The most successful attempt of 
the application of fractals to image compression has been 
by Iterated Systems Inc. While their methodology is a 
well-guarded trade secret, the authors have tested the re- 
sults of their fractal compression scheme against the JPEG 
baseline algorithm. This will be discussed in more detail 
in Section IV. 


III. JPEG Algorithms 

With the advent of multi-media services offered by the 
64-Kbit /sec Integrated Services Digital Networks (ISDN), 
there is a strong urge to define a standard for applica- 
tions as diverse as photo- videotex, desktop publishing, 
graphic arts, color facsimile, photojournalism, medical sys- 
tems, and many others. The JPEG was formed under 
the joint auspices of the International Standards Organi- 
zation (ISO) and the Comite Consultatif International de 
Telephone et Telecommunication (CCITT) at the end of 
1986 for the purpose of developing an international stan- 
dard for the compression and decompression of continuous- 
tone, still-frame, monochrome, and color images. 

The JPEG-proposed algorithm has three major com- 
ponents. The first is a baseline system that provides a 
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simple and efficient algorithm that is adequate for most 
image-coding applications. The second is a set of extended 
system features that allows the baseline system to satisfy a 
broader range of applications. Among these optional fea- 
tures are 12-bit/pixel input, progressive sequential and hi- 
erarchical build-up, and arithmetic coding. The third is an 
independent, differential, pulse-code modulation (DPCM) 
scheme for applications that require lossless compression. 
The following is a brief description of the JPEG baseline 
system only. The reader is referred to the JPEG proposal 
[4] for a complete description of the components and the 
algorithms. 

The JPEG baseline system is a transform-based algo- 
rithm consisting of three stages. The first stage is a dis- 
crete cosine transform (DCT). The output of the DCT is 
then quantized and in the final stage the quantized output 
is encoded by variable length codes. The original image is 
partitioned into 8x8 pixel blocks and each block is inde- 
pendently transformed by the DCT. The transform coeffi- 
cients are then quantized using a user-defined quantization 
template that is fixed for all blocks. Each component of 
the quantization template is an 8-bit integer and is passed 
to the receiver as part of the header information that is re- 
quired for every image. Up to four different quantization 
templates can be specified; for example, different quan- 
tization templates may be used for the different compo- 
nents of a color image. The JPEG baseline system supplies 
two default quantization templates: one for the luminance 
component (the T-component) and the other for the two 
chrominance components (the 7 and Q components). The 
top-left coefficient in the two-dimensional DCT block [i.e., 
the (0, 0) coefficient] is the DC coefficient and is propor- 
tional to the average brightness of the spatial block. The 
remaining coefficients are called the AC coefficients. After 
quantization, the DC coefficient is encoded with a lossless 
DPCM scheme using the quantized DC coefficient from 
the previous block as a one-dimensional predictor. For the 
baseline system, up to two separate Huffman tables for 
encoding the resulting differential signal can be specified 
in the header information. A default Huffman table for 
DC encoding is given in the JPEG proposal. The encod- 
ing of the quantized AC coefficients uses a combination 
of runlength and Huffman coding techniques. There are 
many zeros in the quantized AC coefficients, especially in 
the high frequencies. The AC coefficients that are close 
(respectively, far) in location to (0, 0) are the low (respec- 
tively, high) frequencies. Typically, high frequencies have 
low energies. The two-dimensional block of quantized coef- 
ficients is transformed into a one- dimensional vector using 
a zigzag reordering so that the coefficients are arranged 
in approximately decreasing order of their average energy. 
This creates a combination of nonzero values at the begin- 


ning of the vector and long runs of zeros thereafter. To 
encode the AC coefficients, each nonzero coefficient is first 
described by a composite 8-bit value, denoted by 7, of the 
form (in binary notation) 

7 = NNNNSSSS 

The four least significant bits, SSSS } define a category 
for the coefficient amplitude. The values in category k 
are in the range (2*- 1 ,2* - 1) or (-2 k + l 5 -2 fc ~ 1 ), where 
1 < Ar < 10. Given a category, it is necessary to send 
an additional k bits to completely specify the sign and 
magnitude of a coefficient within that category. The four 
most significant bits in the composite value, i.e., N N N N , 
give the position of the current coefficient relative to the 
previous nonzero coefficients. The runlengths specified by 
NNNN can range from 0 to 15, and a separate symbol 
7 = 11110000 = 240 is defined to represent a runlength 
of 16 zeros. In addition, a special symbol, 7 = 0, is used 
to code the end of a block (EOB), which signals that all 
the remaining coefficients in the block are zero. Therefore, 
the total symbol set contains 162 members (10 categories 
x 16 runlength values -f- 2 additional symbols). The out- 
put symbols for each block are then Huffman coded and 
are followed by the additional bits required to specify the 
sign and exact magnitude of the coefficient in each of the 
categories. Up to two separate Huffman tables for the AC 
coefficients can be specified in the baseline system. A de- 
fault runlength/Huffman table for AC encoding is given in 
the JPEG proposal. 


IV. The Comparison 

For the comparison of the fractal and the JPEG algo- 
rithms, a set of ten 320 x 200 8-bit grey-scale images in 
Sun raster format was sent to Iterated Systems Inc. Some 
of these images are often used for image-compression ex- 
periments, and others are planetary images. The complete 
list is (1) an air scene, (2) baboons, (3) a couple in a room, 
(4) Lena (portrait of a young woman often used in image 
compression tests), (5) peppers, (6) light effects pattern, 
(7) the moon, (8) Miranda, and (9) and (10) two images 
of the planet Saturn. Each of the images was compressed 
to eleven files of sizes ranging from 3 to 13.5 Kbytes, 
The reconstruction was done by the decompression soft- 
ware P.OEtf™ Developers 3 Kit— Grayscale Still devel- 
oped by Iterated Systems Inc. for this purpose and the 
format conversion developed by one of the authors. The 
reconstructed images were then compared to the original 
to obtain the mean-square-error (MSE) and the peak-SNR 
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values versus the number of bits per pixel (i.e., compres 
sion ratio). Peak SNR is defined as 101og 10 (255 f/MSE, 
and is measured in dB’s. Similar curves were obtained for 
the JPEG baseline algorithm as implemented at JPL [5]. 
The results are given in Figs. 1 and 2 for images (5) and 
(8). It is clear from the figures that the JPEG algorithm is 
superior to the fractal algorithm by as much as 6 dB at low 
compression ratios and 2 dB at high compression ratios. 
The curves presented in this article contain the range of 
compression ratios obtained by Iterated Systems Inc. for 
the ten images referred to earlier. It should be pointed out 
that no conclusion regarding the effectiveness of the fractal 
method at higher compression ratios is warranted at this 
time. Qualitatively, at high compression ratios of about 
16:1 to 20:1, the fractal scheme exhibits solar flare, while 
the JPEG baseline algorithm suffers from the tiling effect. 


No post-processing was done with the application of the 
JPEG algorithm. The authors are unaware of any post- 
processing in the fractal algorithm. The tiling effect can 
be somewhat alleviated by the application of certain filters 
which involve weighted averages of nearby neighbors. The 
authors know of no method for dealing with the solar flare 
effect. 


V. Conclusion 

A comparative study of the JPEG baseline algorithm 
and the fractal method was conducted by the authors. In 
terms of the mean square error and peak SNR, the JPEG 
algorithm was superior in every case. 
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* FRACTAL COMPRESSION 
° JPEG COMPRESSION 
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Fig. 1. Fractal performance on peppers. 
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Fig. 2 . Fractal performance on Miranda. 



